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ABSTRACT 

We have performed accurate iron abundance measurements for 44 red giants (RGs) in the Carina 
dwarf spheroidal (dSph) galaxy. We used archival, high-resolution spectra (R ~ 38 000) collected with 
UVES at ESO/VLT either in slit mode (5) or in fiber mode (39, FLAMES /GIRAFFE- UVES) . The 
sample is more than a factor of four larger than any previous spectroscopic investigation of stars in 
dSphs based on high-resolution (R > 38 000) spectra. We did not impose the ionization equilibrium 
between neutral and singly-ionized iron lines. The effective temperatures and the surface gravities 
were estimated by fitting stellar isochrones in the V, B-V color-magnitude diagram. To measure 
the iron abundance of individual lines we applied the LTE spectrum synthesis fitting method using 
marcs model atmospheres of appropriate metallicity. For the 27 stars for which we measured both 
Fei and Fen abundances, we found evidence of NLTE effects between neutral and singly- ionized iron 
abundances. The difference is on average ^0.1 dex, but steadily increases when moving from the 
metal-rich to the metal-poor regime. Moreover, the two metallicity distributions differ at the 97% 
confidence level. Assuming that the Fen abundances are minimally affected by NLTE effects, wc 
corrected the Fe I stellar abundances using a linear fit between Fe I and Fe n stellar abundance deter- 
minations. We found that the Carina metallicity distribution based on the corrected Fei abundances 
(44 RGs) has a weighted mean metallicity of [Fe/H] = —1.80 and a weighted standard deviation of 
a = 0.24 dex. The Carina metallicity distribution based on the Fen abundances (27 RGs) gives sim- 
ilar estimates ([Fe/H] = —1.72, a = 0.24 dex). The current weighted mean metallicities are slightly 
more metal poor when compared with similar estimates available in the literature. Furthermore, if 
we restrict our analysis to stars with the most accurate iron abundances, ~ 20 Fei and at least three 
Fen measurements (15 stars), we found that the range in iron abundances covered by Carina RGs 
(~ 1 dex) agrees quite well with similar estimates based on high-resolution spectra. However, it is 
a factor of two/three smaller than abundance estimates based on the near- infrared Calcium triplet. 
This finding supports previous estimates based on photometric metallicity indicators. 
Subject headings: galaxies: dwarf — galaxies: individual (Carina) — galaxies: stellar content — stars: 
abundances — stars: fundamental parameters 
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1. INTRODUCTION 

The Carina dSph galaxy is a fundamental benchmark 
to constrain the formation and evolution of dwarf galax- 
ies and a key laboratory to improve current knowl- 
edge of low- and intermediat e-mass stars' evolutionary 
and pulsation properties (e.g. iSmecker-Hane et a,l.lll999t 
iMonelli et aLll2003t iBono et al.ll2010D . Quantitative con- 
straints on its stellar content require not only precise 
multiband photomet ry over the entire body of the galaxy 
(|Stetson et al.H20lU) , but also accurate measurements of 
the mean metallicity and of the spread, if any, in metal- 
licity. 

Using the multi-fiber sp ectrograph ARGUS at the 
CTIO 4 m Blanco telescope, ISmecker-Hane et all (|1999l ) 
collected low-resolution (R ~ 3 000) spectra for 52 Red 
Giants (RGs), covering the near-infrared (NIR) Cal- 
cium triplet (CaT) region, and found a mean metal- 
licity of [Fe/H] = —1.99 and a small spread in metal- 

licity (cr[Fe /Hl = 0-25 dex). On the other hand, 

iKoch et al.l ((2006) collected medium-resolution spectra 
(R ~ 6 500) for 437 RGs wit h the multi-fiber sp ec- 
trograph FLAMES/GIRAFFE (|Pasquini et all 12002 ) in 
MEDUSA mode at ESO/VLT,and using the same diag- 
nostic found a mean metallicity of [Fe/H ] = —1.90 ±0.01 
(metallicity scale bv lZinn fe~W cst 1984|) and a significant 
spread in metallic ity (full range —3.0 < [Fe/H] < 0.0). 
IKoch et al.l (|2006[ ) transformed the individual equivalent 
widths (EWs) of the CaT line s into iron abundances us- 
ing the calibration provided bvlRutledee et al.l (119971) . In 
a more recent investigation, HclmLetalJ (|2006|) studied 
the same spectroscopic data and found a similar mean 
metallicity, ([Fe/H] = —1.7 ±0.1), but a smaller spread 
in iron abundance (—2.3 < [Fe/H] < —1.3). However, 
they used a different c alibration between t he EW of CaT 
lines and metallicity (jTolstov et al.l [200lT) and different 
criteria for selecting candidate Carina stars. The differ- 
ence between the two metallicity distributions appears as 
a steeper metal-rich tail in the latter distribution when 
compared with the former one. 

High-resolution spectra (R ~ 40 000) are available for 
a sample of five bright RG stars, collected with the slit 
mode of Ultraviolet a nd Visual Echelle Spectrograph 
(UVES) at ESO/VL T dDekker et al l lpOOOl) . The mean 
metallicity found by iShetrone et al.l (|2003l) from these 
spectra is [Fe/H] = —1.64 with crpe/H] = 0.4 dex. The 
spread decreases to 0.25 dex if we remove the most metal- 
rich star in the sa mple. Independent measurements by 
IKoch et alJ (|2008f ). using high-resolution (R ~ 38 000) 
spectra for ten bright RGs collected with the red arm of 
UVES from FLAMES/GIRAFFE fibers, give a similar 
mean metallicity of [Fe/H] = —1.69, but a spread that is 
a factor of two larger (crpe/H] = 0.51 dex). 

A detailed spectro scopic analysis of Car ina RGs was 
recently provided by iLemasle et al.1 (|2012f) . They used 
FLAMES/GIRAFFE spectra for 35 RGs collected with 
grisms HR10 (R ~ 19 800), HR13 (R ~ 22 500) and 
HR14 now (R ~ 17 700). They found a mean abun- 
dance based on Fei lines (^[Fei/H] = —1-63 ± 0.01, 
o'lFoi/H] — 0.27 dex, weighted mean) that agrees quite 
well with similar estimates available in the literature. On 
the other hand, the mean abundance based on Fen lines 
(M[Feii/H] = -1.24 ±0.01, fT[Fcii/H] = 0.31 dex, weighted 
mean) is slightly more metal-rich, but the difference is of 



the order of la. Moreover and even more importantly, 
they found that their metallicity estimates range from 
[Fei/H] = -1.2 to [Fei/H] = -2.5, thus supporting a 
corresponding decrease in the metallicity spread. 

A more recen t analysis has been provided by 
IVenn et al.l (I20T1 ). using high-resolution spectra col- 
lected either with FLAMES / GIRAFFE- UVES (seven 
RGs) or with Magellan-MIKE (R ~ 28 000, blue, ~ 
22 000, red, two RGs). They found, by using EWs, 
mean abundances based on Fei (^[Fci/H] = —1-79 ±0.09, 
o'pci/H] — 0.54 dex, weighted mean) and on Fen lines 
(/^[Fcii/Hl = —1-76 ±0.09, <T[Feii/H] = 0-53 dex) that agree 
quite well with similar estimates available in literature. 

Current spectroscopic investigations support the ev- 
idence that Cari na may host very metal-poor stars 
([Fe/H] < -3.0, iHelmi et al.l [20061) . This evidence, 
if supported by independent spectroscopic investiga- 
tions, would imply that Carina underwent a fast chem- 
ical enrichment starting in a metal-poor environment 
and rapidly approaching a metal-intermediate regime 
([Fe/H] ~ -1.6). 

Photometric investigations agree with the spectro- 
scopic measurements concerning the mean metallicity, 
but suggest a small sp read in metallicity (crpc/H] ~ 
0.25 dex, IBono et al.l [20TrH 

°"[Fe/H] ~ 0.35 dex, 
ILianou et al.ll20rT[ ). In this work we per form a reanaly- 
sis of th e Carina sp e ctra in vestig ated by She trone et al.l 
( 2003), IKoch et al.l (l200l and IVenn et al.1 (|2012h to- 
gether with the analysis of a new high-resolution spec- 
troscopic data set. In total we have 89 spectra for 72 
stars. However, 18 of them are probable non-members, 
two of them are carbon stars, for six of them we have not 
been able to measure the radial velocity and for two of 
them we have not been able to measure the iron abun- 
dances. Overall, we provide new homogeneous iron abun- 
dance measurements for 44 Carina RGs based on high- 
resolution spectra collected with UVES (slit mode, five) 
and with FLAMES / GIRAFFE-UVES (multifiber mode, 
39) at ESO/VLT. 

This is the largest sample of high-resolution (R > 
38 000) spectra ever collected for a dSph galaxy. A 
similar spect r oscop ic approach was also adopted by 
iLetarte et all (|2010D for RGs in the Fornax dSph. They 
used FLAMES/GIRAFFE spectra for 81 RGs collected 
with grisms HR10, HR13 and HR14 ncw . A set of spec- 
tra were also collected with the old version of the grism 
HR14 id (R ~ 28 800), but they were rescaled to the res- 
olution of the spectra collected with the grism HR14 new . 

2. OBSERVATIONS AND DATA REDUCTION 

The high-resolution spectra for Carina RGs adopted 
in this investigation were retrieved from the ESO Sci- 
ence Archive. We selected spectra from four differ- 
ent ESO/VLT observing programs collected with either 
UVES (slit mode, nine) or FLAMES / GIRAFFE-UVES 
(multifiber mode, 80). 

The oldest (2000-2001) data seQ includes nine spectra 
of five RGs with visual magnitudes ranging from 17.62 
to 17.92 mag. These spectra were collected with red arm 
of UVE S, include 37 orders an d have already been ana- 
lyzed bv IShetrone et all (|2003l ). The second (2003) data 

18 ESO programs 065.N-0378( A) 066.B-0320 ( A), FT : Tolstoy. 
These are the spectra adopted by Shetrone et al. (2003). 
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selF^I comprises individual spectra for 33 RGs collected 
with the FLAMES /GIRAFFE-UVES red arm; their vi- 
sual magnitudes range from 17.32 to 18.98 mag. The 
orders and the wavelength coverage of these spectra are 
the same as the UVES spectra. They have already been 
analyzed by IKoch et al.l ()2008D who found that two of 
them are candidate carbon stars (stars Car54, Car55, 
or using Kock's IDs, LG04a_000057 and LG0 4b.000569). 
We su pport the classification suggested by IKoch et all 
(2008), and have not included these stars in the current 
abundance analyses. The third (2005) data setP 7 ! and 
collected with FLAMES/GIRAFFE- UVES. The visual 
magnitudes of the RGs range from 17.65 to 18.00 mag. In 
the fourth (2007-2008) data seQ there are 40 spectra for 
40 different stars collected with FLAMES/GIRAFFE- 
UVES, and their visual magnitudes range from 17.61 to 
18.68 mag. 

All the spectra collected with the red arm of UVES 
cover the wavelength range 4780 - 6825 A and the 
signal-to-noise ratio (S/N) is typically of the order 
of 30 (A - 6750 A). The spectra collected with the 
FLAMES /GIRAFFE-UVES red arm cover the same 
wavelengths and have S/N (at A ~ 6750 A) ranging from 
~15±5 (21 stars) to ~30±10 (11 stars) and -45±5 (7 
stars). 

We ended up with a sample of 89 high-resolution spec- 
tra for 72 stars, located across Carina's central regions 
and covering the bright portion of the RG branch (see 
Fig. 1). We chose to use spectra taken with the red arm 
of UVES (centered at 5800 A) because it is minimally 
affected by contamination from sky lines. The UVES 
spectra were reduced using IRApQ and extracted using 
the IRAF task apall. Wavelengt h calibration was pe r- 
formed using reference lines from (|Murphv et al.ll2~007l ). 

3. RADIAL VELOCITY AND PHOTOMETRIC ANALYSIS 

We measured the radial velocity (RV) of each spec- 
trum followin g the same approach as adopted by 
iFabrizio et al.l (|2011l) . using two dozen heavy-element 
lines ranging from 6136 to 6200 A. We measured the RV 
of 66 stars (81 spectra). For six stars the measurement 
of the RV was not possible due to the poor quality of 
the spectra. Among the stars with RV measurements, 
46 (61 spectra) are candidate Carina stars (with 212 < 
RV < 243 km s" 1 . IFabrizio et al.l 1201 ll ). The others are 
either candidate field stars (18, RV < 100 km s _1 ) or 
candidate carbon stars (two) and they were not included 
in the current analysis. 

To provide accurate estimates for the stellar param- 
eters of the spectroscopic targets (effective tempera- 
ture T c g and surface gravity lo gq), we use d the multi- 
ba nd photometry discu ssed in IBono et"aT1 (|2010l ) and 
in iStetson et all (|2011f ). Moreover, we adopted dif- 
ferent scaled-Solar cluster isochrones from the BaSTI 

19 ESO programs 17 1.B-0520(A) ( B)(C) , PI: Gilmore. These are 
the spectra adopted by Koch et al. (2008). 

20 ESO programs 074.B-0415(A), 076.B-0146(A), PI: Tolstoy. 
These are the hig h-resolution spectra for seven RGs adopted by 
IVenn et al.l ( 12013) 

21 ESO programs 180.B-0806(B), PI: Gilmore 

22 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



data baseH (|Pietrinferni et all 12004 120061) . We selected 
isochrones at fixed age (12 Gyr) for three different chem- 
ical compositions ([Fe/H] = -1.50,-1.79,-2.27). We 
also adopted a true distance modulus D M n = 20.15 mag 
and a reddening E(B — V) = 0.04 mag dDaH'Ora et al.1 
2003). Finally, we performed a linear regression among 
the different isochrones and derived analytical relations 
connecting the visual magnitude V, the B-V color and 
either the surface gravity or the effective temperature. 
On the basis of the above relations and of the observed 
magnitudes and colors, we estimated the atmospheric pa- 
rameters for each spectroscopic target. The adopted ap- 
proach and the analytical relations will be described in 
a future paper (Ferraro et al. in preparation). 

The top panel of Fig. 1 shows the position of the 
spectroscopic targets in the V, B-V Color-Magnitude 
Diagram (CMD) to gether with the entire photometric 
catalog (grey dots, IBono et all 120101 ). The black dots 
mark the 46 candidate Carina stars, while the tiny or- 
ange dots display the 20 candidate field stars and the 
two large orange dots show the carbon stars. The blue 
sq uares and the red s tars s how the high-re s olutio n sample 
of IKoch et al.l (|2008l ) and iShetrone et ahl (|2003t ). respec- 
tively. The solid lines display the adopted isochrones. 
Taking into account current uncertainties in the Ca- 
rina true d istance modulus, the reddening, the mean 
metallicity (iDall'Ora et alj|2003t iPietrzvnski et al.ll2009t 
IBono et al.ll201Ct iLemasle et al.ll2012h and the spread in 
age jstetson et al.ll2011l ). we performed a series of simu- 
lations and found that the typical uncertainties in tem- 
perature and gravity are £T off ~ 70 K and ei ogg ~ 0.2 dex 
(see Table 1). The bottom panel of Fig. 1 shows the same 
isochrones as the top panel, but in the log g vs T e g plane 
together with the target stars. The error bars plotted 
in the bottom right corner represent the aforementioned 
uncertainties. 

4. SPECTROSCOPIC ANALYSIS 

The iron abundance analysis was performed following 
the classical spectrum-synthesis method for both Fe i and 
Fen lines, but with one difference: we did not impose 
LTE ioniz ation equilibrium b etween the Fei and Fen 
lines (e.g. iKraft fc Ivaf3l200l . This means that we trust 
the surface gravity determined from the optical photom- 
etry and cluster isochrones. For the highest S/N spectra 
collected with UVES, we selected Fe i and Fen lines from 
the valiQ data base (jKupka et al.l l2"000). We ended up 
with a list of 123 Fei and 18 Fen lines in the wavelength 
range covered by our spectra. Among them, 45 Fe i and 
5 Fen lines are located in overlapping orders. 

By using the photometric estimates of both T e g and 
logff, and the Ca rina mean metallicity [Fe/H] = —1.70 
(jKoch et all 120081). we interpolated t he marcS model 
atmospheres (fGustafss onet al.l 120081 ) with a modified 
version of the Masseron (2006) codcQ for each Carina 
star. The individual models assumed spherical geome- 
try, an a-enhanced (Ja/Fe] = +0.4) chemical mixture, 
a mass value of lMgO and a constant microturbulence 

23 Available at the URL : Ihttp : / /albione .oa-teramo. inaf . it / 1 

24 Available at the URL: http://vald.astro.univie.ac.at 

25 Available at the URL: http://marcs.astro.uu.se 

26 Available at the URL: http://marcs.astro.uu.se/software.php 

27 Note that for stellar masses ranging from 0.8 to 1.2Mq, the 
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velocity (£ = 2 km s _1 ). The value adopted for the mi- 
croturbulence velocity follows the estimates provided by 
iShetrone et al.l (|2003( ) for Carina red giants using metal- 
poor stellar model atmospher es . Simila r estim ates for 
the £ values were provided by iTheveninl (|1998f ) for RG 
stars with similar atmospheric parameters in the globular 
uj Cen. The Carina targets span a modest range in at- 
mospheric parameters, therefore we decided to adopt the 
same £ value for the entire sample. The synthetic spectra 
were computed using the ID, plane-parall el, LTE, radia - 
tive transfer code moogP^I (2009 version, ISnededHom 
Note that plane-parallel radiative transfer for giant stars 
is more accurate when us ing spherical model atmospheres 
(|Heiter fe Eriksson|[2006l ). 

We chose the So lar chemical composition from 
iGrevesse et al.1 (|2007f ). i.e., a Solar iron abundance of 
Ap e = 7.4E0 Oscillator strengths / for Fei and 
Fen lines are from vald. When available, the pa- 
rameters for collisional dampin g with neutral hydrogen 
are based on quantum theory (lAnstee fe Q'Maral[T995t 
iBarklem fe O'Maral 119971 : IBarklem et al l 11998ft : other- 
wise the Unsold classical recipe (lUnsoldlll955h~ was used 
with an enhanced damping factor of Fjj = 1.5. We note 
that calculations based on quantum theory are available 
for all the selected Fen lines. The strategy adopted for 
the iron abundance analysis was dictated by the quality 
of the observations: high-resolution (R ~ 38 000) and a 
reasonable signal-to-noise ratio (15 < S/N < 45). We 
decided to use an eye fitting procedure rather than an 
automatic EW fitting procedure, since most of the weak 
lines are either noisy or blended. 

The selected Fe i and Fe n lines were normalized using 
a linear fit of the local continuum identified by eye in- 
spection. Using the moog code, we synthesized spectra 
with different metallicities. The synthetic spectra were 
convolved with a Gaussian broadening function to repro- 
duce the instrumental resolution. We neglected the effect 
of stellar rotation. The fit of the lines (45 Fe i + 5 Fe n) 
that appear in two overlapping orders were treated in- 
dependently and the iron abundance for these lines was 
estimated as a weighted mean. After a preliminary anal- 
ysis of the stars with multiple spectra, we decided to ne- 
glect the spectra with lower S/N and we only analyzed 
the best available spectrum for each star (46 spectra). 

It is worth mentioning that in the abundance analysis 
we did not include lines with A < 4920 A to overcome 
subtle uncertainties with the coherent scattering by hy- 
drogen in the contin uum that was not included in moog 
(|Sobeck et al.l 1201 lh . To constrain quantitatively the 
possible error in abundance in this wavelength range we 
computed the abundances of the three Fei lines at 4924, 
4950 and 4973 A. We used a marcs model atmosphere as- 
suming T e ff = 4500 K, log g = 1 d ex and [Fe/H] = -1.5 
with the code TURBOSPECTRUivPlf Alvare z fe Plezlll998l ) 
which accounts for the coherent scattering by hydrogen in 
the continuum. We found no significant difference with 
the abundance measurement based on the moog code, 
and indeed the differences in the cores of the lines repre- 



sent less than five percent of the line equivalent widths. 

We performed a first preliminary mean iron abun- 
dance analysis for each star. Subsequently, we com- 
puted a new set of model atmospheres using the same 
values of T c g and log g, but using the new individ- 
ual iron abundances. We also decided to use scaled- 
Solar Qa/Fe] = 0.0) marcs models, since detailed abun- 
dance measurements of RGs in nearby dwarf galaxies 
(jShetrone et alj 120031: Koch et al.l 120081 ) indicates that 
they show only a mild enhancement in a-elemcnts. For 
stars with [Fe/H] < —2.0, we used a-enhanced models 
because the very metal-poor, scaled-Solar marcs models 
are not available. In this new iron abundance analysis, 
we double checked all the iron measurements and rejected 
strong, noisy, blended or too weak Fe i and Fe n lines from 
the analysis. The adopted line list with wavelengths for 
both neutral and ionized iron lines is presented in Ta- 
ble 2 together with ion identification, excitation poten- 
tial (EP) and log gf values. We ended up with 83 Fei 
and 11 Fen lines. Among them, 25 Fei and 4 Fen lines 
are located in overlapping orders. We performed several 
experiments with synthetic spectra and we found that 
the accuracy of the line fits ranges from 0.1 to 0.15 dex. 

The atmospheric uncertainty in the line-by-line abun- 
dance measurements was estimated usi ng the formul a 
based on the EW uncertainty suggested bv lCavrel (|1988j) . 
In particular, we adopted a mean equivalent width EW 
(typical values of the selected lines are ~ 50 mA for Fe i 
and ~ 30 mA for Fen) and the S/N of each individual 
spectral order (estimated using the IRAF task splot). 
The final iron abundance for each star was computed as 
the weighted mean of the line-by-line abundances. Indi- 
vidual weights were chosen as the inverse of the atmo- 
spheric uncertanties to limit the influence of the outliers. 

For each star, in Table 2 are also listed the line-by- 
line EWs. The EWs were measured using the IRAF task 
splot, assuming a single Gaussian line profile, since the 
selected lines are not blended nor asymmetric. Typi- 
cal uncertainty on the widths is of the order of 5%, 
and it is caused by the estimate of the continuum. It 
is worth mentioning that the EWs listed in Table 2 
are only given for completeness of the current spectro- 
scopic analysis. The iron abundances we provide are 
only based on the spectrum synthesis fitting method. 
The differences betwe en our EWs and those m easured by 
IShetrone etHI 1)20031 ) and IKoch et al.l (|2008D are min- 
imal (weighted mean, EW our — EW aut hor ~ —5.4 and 
~ 1.7 m A, respectively) . The difference with EWs pro- 
vided by iLemasle et al.l (|2012l ) is larger (eleven stars in 
common, ~ -22 mA). The comparison between the 
EWs provi d ed by ILemasle et al.l ^2012) and those by 
IKoch et all ([2008) gives a similar difference (one star in 
co mmon, ~ —27 mA ) . The difference with EWs provided 
by I Venn et al.l (|2012t ) is quite smaller (seven stars in com- 
mon, ~ —13 mA) , and compared with EWs provided by 
IKoch et all (2008) gives a similar difference (three stars 
in common, ~ — 18 mA). 



change in logg, at fixed T eff , is smaller than 0.1 dex 5. RESULTS AND DISCUSSION 

28 Available at .the URL:[^://ww as.utexas.edu/~chris/moog.html| 5 ^ Detailed Fe I and Fe II stellar abundances 
A-p c = log(Np c /Nu) + 12 where Np e and A<h are the number 

densities of Fe and H We measured [Fei/H] for 44 stars and [Fen/H] for 27 

30 Available at the URL: |htt P ://www.graal.univ-mont P 2.fr/hosted/plei| t- iairR ou t f the 46 stars in our sample. On average, 
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the former are based on more than twelve measurements 
(five stars have a single measurement), and the latter 
are based on more than three measurements (eight stars 
have a single measurement). The iron abundance was 
not measured in two metal-poor stars (Carl7, Car26) be- 
cause the lines available were noisy, blended or misiden- 
tified. The individual [Fei/H] and [Fen/H] abundance 
measurements and their atmospheric errors for the en- 
tire sample are listed in Table 2. The mean [Fei/H] and 
[Fcii/H] abundances for the entire sample, with their 
weighted standard deviations are given in columns 10 
and 12 of Table Columns 11 and 13 list the num- 
ber of Fe i and Fe n lines used to estimate the individual 
weighted mean abundances. 

We carried out a more detailed analysis to constrain 
the impact that the adopted atmospheric parameters 
have on individual abundance measurements. To address 
this classical issue, we computed synthetic spectra with 
a series of model atmospheres at fixed chemical compo- 
sition, changing the effective temperature by 100 K, the 
surface gravity by 0.3 dex and the microturbulcnce ve- 
locity by 0.5 km s" 1 . Note that we are using generous 
estimates for the uncertainties affecting the atmospheric 
parameters (see §3). In particular, we selected the star 
Car23, since its effective temperature (~ 4400 K) and 
surface gravity (~ 0.80 dex) can be considered as average 
values for the entire sample. The results listed in Table 3 
indicate that a change of 100 K in effective temperature 
causes a change in the mean Fei and Fen abundances 
of the order of 0.1 dex. On the other hand, a change of 
0.3 dex in logg has a minimal impact on the Fei abun- 
dance and a difference of the order of 0.15 dex for the 
Fen abundance. The differences in abundance caused 
by a change of 0.5 km s _1 in microturbulence velocity 
are of the order of 0.05 and 0.07 dex for the Fei and 
Fe ii abundances, respecti vely. These findin gs support 
the analysis performed by Ko ch et al.l ([2008) and indi- 
cate that plausible changes in the atmospheric parame- 
ters introduce changes in the inferred iron abundances 
smaller than 0.15 dex. To account for the uncertainties 
affecting the atmospheric parameters, the errors on the 
weighted mean metallicities listed in column 10 and 12 
of Table 1 were estimated by summing in quadrature the 
weighted standard deviation with the three typical errors 
stated above. 

The iron abundances of the 27 stars for which we have 
both [Fei/H] and [Fen/H] measurements (Table 1) in- 
dicate that the [Fei/H] abundances are systematically 
lower than [Fen/H]. The mean difference ranges from 
~0.1 to 0.2 dex (Carl2, Car3). Five stars show an op- 
posite trend, but their Fe n measurements are based on 
one or two lines, and two out of the five have Fe i mea- 
surements based on two lines. 

The detailed iron abundance analysis for five Carina 
stars, covering the entire range of apparent magnitudes, 
is shown in Fig. 2. To validate the accuracy of our T e s 
determinations, we plotted for each star the abundance 
of individual lines as a function of the Excitation Po- 
tential (EP) of the lines. We performed linear fits and 
found that the slopes are within ±0.06 dex/eV (see la- 
beled values in Fig. 2), except for a few stars with a 

31 Columns 10 and 12 give the atmospheric uncertainty for the 
stars with only one line measurement. 



limited number of Fe i measurements. We tested the Fe i 
lines (left panels) , since they are based on a larger num- 
ber of measurements. The data plotted in Fig. 2 indicate 
that weighted standard deviations for Fe i and Fe n lines 
are on average 0.12 and 0.18 dex, respectively. We also 
checked the impact of a change in microturbulence ve- 
locity on star Carl4 by assuming £ = 2.5 km s . The 
slope changed by less than 0.02 dex/eV, since the satu- 
rated Fe i lines were not taken into account in the current 
analysis. 

5.2. Comparison with previous investigations 

To validate our abundance analysis, we performed a 
detailed comparison with similar measurements available 
in the literature. We defined the iron abundance differ- 
ence A[Fe/H] = [Fe/H] our -[Fe/H] author after shifting the 
[Fe/H] author to the Solar iron abundance that we used 
(Af e = 7.45. iGrevesse et al.l 120071 ). The data plotted in 
Fig. 3 show A[Fei/H] (left) and A[Fen/H] (right) versus 
our iron abundance. The vertical error bars are the sum 
in quadrature of our atmospheric dispersion and the er- 
rors given by the different authors, while the dotted lines 
display the 1 a difference. 

We have five stars in common with She trone et al.l 
(p003T) (Car2, Car3, Car4, CarlO, C ar 12). Their Solar 
iron a bundance reference is 7.52 from IGrevesse fe Sauvall 
(1998). The data plotted in the panels a) of Fig. 3 
show that our iron abund ances are systematicall y more 
metal-poor than those by iShetrone et al.l (|2003D . espe- 
cially for Fei. The weighted mean and weighted stan- 
dard deviation of the difference are Ma[FcI/h] — —0.25 ± 
0.07, cta[FoI/h] = 0.17 and /UA[Feii/Hi = -0.17 ± 0.11, 
°~A [Feii/H] = 0.14 dex. Star Car SEE the object with 
the largest discrepancy, suggests that the iron abun- 
dance from these author s is an overestimat e. This hy- 
pothesis is supported by iKoch et ail (|2006[ ) who found 
[Fe/H] = —1.93 after accounting for the difference in 
the zero-point. This agrees better with our abundances 
([Fei/H] = -2.14, [Fen/H] = -1.97) than with the 
Shetrone's estimates ([Fei/H] = -1.58, [Fen/H] = 
— 1.56, using our iron Solar abundance). Neglecting 
this star, our abun dances agree quite well with the 
IShetrone et alj (|2003t ) measurements, and indeed the dif- 
ferences decrease to /^A[Fei/H] = — 0.17±0.08, cta[foI/h] = 
0.07 and HA[FeU/ti\ = -0.09±0.12, ct A [fcII/h] = 0.04 dex. 
We adopted the same moog radiative transfer code and 
the same marcs model atmospheres (values of T e g and 
log<? similar), but we used more recent versions (2009 for 
moog and 2008 for marcs). The micr oturbulence veloc- 
ities adopted bv IShetrone et al.1 (|2003l ) range from 1.9 to 
2.2 km s^ 1 and are quite similar to the value we adopted. 
We also performed a series o f tests to constrain the dif- 
ference between the results of lShetrone et al.l (|2003f l and 
the current abundances. The differences of EWs mea- 
sured for this work (s ee Table 2) and those published by 
IShetrone et alj (|2003l see their Table 4) is minimal (a 
few mA, see §4). The oscillator strengths the y adopted 
come from the papers of the Lick- Texas group (Fulbright 
120001 and references therein) and from the NiSTon-line 

32 The reader interested in more detailed analysis of the pho- 
tometric and spectroscopic properties of Car3 is referred to the 
Appendix. 
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Atomic Spectra DatabaseF^I. We have already mentioned 
that our atomic data come from the vald database, and 
it is very difficult to critically evaluate the accuracy of 
atomic data coming from different databases. This is 
the reason why we give a very low statistical weight to 
the stars in our sample with iron abundances based on 
only one or two lines. This evidence indicates that the 
discrepancies might result from the different approaches 
adopted to estimate the abundances. Their abundances 
were estimated using the EW fitting method and forc- 
ing ionization equilibrium between Fe i and Fe u, whereas 
we adopted spectrum synthesis fitting method and did 
not force ionization equilibrium. It is also worth not- 
ing that the n u mber of Fe i and Fe n lines adopted by 
IShetrone et al.l (|2003ft is three times larger than our se- 
lection. 

We also have ten stars in common with iKoch et al.1 

(2008); In their spectroscopic analysis (see also 

IKoch et ahl 120061). these authors adop ted the metallic- 
ity scale of iCarretta fe Grattonl (|1997ft . Therefore, we 
assumed that they adopted a Solar iron abundance of 
7.52. The comparison, after the shift in the zero-point, 
is shown in panels b) of Fig. 3. We found that we agree, 
within lcr, with these authors for the Fe n abundance 
with /xa[fcII/h] = —0.21 ±0.08 and cta[FcII/h] = 0-27 dex. 
For Fei the discrepancy is larger, with /iA[Fei/H] = 
—0.34 ± 0.07 and cta [Fei/H] = 0.26 dex. In this case the 
difference in the abundances might be due to a difference 
either in the approach (EW versus spectrum-synthesis 
fitting methods) or in the input phy sics (or both). As 
a matter of fact, IKoch et ail (|2008T ) adopted their os- 
cillator strengths fr om the standard RG star Arcturus. 
Moreover, they used lCastelli fe Kuruczl ((2003) model at- 
mospheres, while we used marcs models and their sur- 
face gravities are on average 0.3-0.5 dex higher than ours. 
The difference in gravity is the consequence of the dif- 
ferent approach in estimating the gravity (forcing the 
balance between Fe i and Fe n versus photometric gravi- 
ties). 

We compared the current Fei and Fen abun- 
dances with s i milar abundances recently provided by 
iLemasle et all (pOll ) from FLAMES /GIRAFFE spec- 
tra with a spectral resolution that is on average a 
factor of two smaller than ours. The data plotted 
in the panels c) of Fig. 3 show that the differences, 
after c orrecting for their refe rence Solar iron abun- 
dance (|Grevesse fe Sauvall [T998). in the Fei abundances 
(a*a[FcI/h] = —0.27 ± 0.09 dex, weighted mean) and in 
the dispersion (<7A[Fei/H] = 0.16) for the eleven stars in 
common are quite similar to the results based on the 
other high-resolution estimates. The difference in the 
Fen abundances is based on only four objects and at- 
tains similar values, namely jUA[Feii/H] = —0.61 ± 0.27, 
°A[Feii/H] — 0-25 dex. Note that our [Fen/H] abundance 
for the star with the largest discrepancy is based on only 
three Fen lines. 

Finally, we compared our Fei an d Fe ii abundances 
with similar abundances provided by I Venn et"all (120121 ) 
using EWs of FLAMES / GIRAFFE- UVES spectra. The 
data plotted in the panels d) of Fig. 3 show that the 
Fe i abundances for the seven stars in common, af- 



ter correcting for the ir reference Solar iron abundance 
(jAsplund et all 12009ft . (ma[f c i/h] = -0.37 ± 0.11 dex, 
weighted mean) and their dispersion (o"A[Fei/Hl — 0.19) 
agree quite well with current estimates. The difference 
in the Fen abundances is based on six objects and at- 
tains similar values, namely /^a[FcII/h] = —0.26 ± 0.13, 
OA[FeII/H] =0.11 dex. 

We also compared our iron abundances to similar es- 
timates by IKoch et al.l ([20061 ). These authors performed 
detailed measurements of CaT EW for 437 RGs and 
transformed the reduced EWs in to iron abun d ances , 
using the metallicity scales from [Zinn fe Westl (|1984l 
ZW84) and lCarretta fe Grattonl ((19971, CG97). We com- 
pared their Fe i abundance estimates with our measure- 
ments of Fe i abundances. We have 25 stars with at least 
two iron line measurements in common with this sam- 
ple and we found that the differences in the ZW84 and 
the CG97 metallicity scales are: Ma[Fo/h] = —0.14 ± 
0.04 (cr A[Fc/H] = 0.15) and AiA[Fc/H] = -0.33 ± 0.04 
("ATFe/Hl = 0.18) dex (see Fig. 4). Our iron abundances 
are slightly more metal-poor, but within 1 a and in better 
agreement with the ZW84 than with the CG97 metallic- 
ity scale. The difference of 0.2 de x between the ZW84 
and the CG97 scale is we ll known (jCarretta fe Grattonl 
[19971 iKraft fe Ivandl2003ft . 

5.3. Carina metallicity distribution 

The data plotted in the top panels of Fig. 5 show the 
weighted metallicity distributions of our sample based 
on Fei (left) and on Fen (right) measurements. To 
constrain the metallicity distr ibution, we assign ed to 
each star a Gaussian kernel (|Di Cecco et al.1 [20101 ) with 
a a equal to the standard deviation of the iron abun- 
dance measurement. The solid lines were computed by 
summing the individual Gaussian over the entire data 
set (see panels a and b in Fig. 5). We found that 
both the median (—1.84 dex) and the weighted mearP^I 
(MFFei/Hl = — 1-90±0.02 dex) of Fe i abundances (44 stars) 
attain similar values. The same outcome applies to the 
Fe ii measurements (27 stars) for which we determined 
M[Fcii/H] = —1-72 ±0.04 dex. The weighted standard de- 
viations attain small values (~ 0.25 dex), but the metal- 
licity distributions cover at least one dex. We also esti- 
mated the metallicity distribution by giving equal weight 
to all stars, and the difference with the weighted mean is 
minimal (< 0.04 dex). 

To constrain the dependence of these results on the 
accuracy of the abundance analysis, we recomputed the 
metallicity distribution using only stars with more than 
twelve Fe i and two Fe n line measurements (see panels c 
and d in Fig. 5). The weighted means and the weighted 
standard deviations reach values similar to those from 
the entire sample, but the range in metallicity decreases 
by at least 0.2 dex. The difference is mainly in the metal- 
poor tail, and indicates that high-quality spectra are re- 
quired to constrain the actual dispersion in iron abun- 
dance of Carina RGs. To further constrain the impact 
of the abundance precision on the metallicity distribu- 
tion we selected the stars with at least three Fen line 



34 The weighted mean metallicitics are estimated by weighting 
the i ron abundance of each star with the quadratic inverse of the 
33 Available at the URL: http://physics.nist.gov/cgi-bin/AtData/mainiBdiilvidual uncertainty. 
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measurements. The Fe i metallicity distribution of the 
same stars is typically based on ~ 20 Fei lines. We 
ended up with a sample of 15 stars, and the weighted 
mean for the Fe i metallicity distribution is slightly more 
metal-rich /xpd/H] = —1-77 ±0.04 dex, while the a [Fei/H] 
decreases to 0.19 dex. The same outcome applies to the 
Fen metallicity distributions. Interestingly enough, the 
Fe i and Fe n metallicity distributions (see panels e and f 
in Fig. 5) cover a range of ~ 1 dex. The use of the best 
data indicates that the spread in metallicity of Carina 
RGs is smaller than previous estimates. Note that this 
finding appears to be minimally affected by selection bi- 
ases, since the metal-poor stars ([Fe i/H] < —2.20) not in- 
cluded in the metallicity distribution of panel e) cover the 
same range in luminosity as the selected starj^l. How- 
ever, these findings call for new abundance analyses to 
constrain the real extent of the metal-poor tail of Carina. 

By taking into account the 27 stars for which we de- 
termined both Fei and Fen abundances, we found that 
the weighted mean based on Fei lines is 0.12 dex lower 
than the weighted mean based on Fen lines. The dif- 
ference brings forward the possible occurrence of NLTE 
effects between Fei and Fen abundances. Our experi- 
ment appears appropriate to constrain this effect, since 
we relaxed the constraint on the ionization equilibrium 
between Fei and Fen lines. This is an interesting find- 
ing, since detailed calculations and the comparison be- 
tween theory and ob s ervations (|Thevenin fe Idiartlll999t 
IThevenin et alJl200H iMashonkina et al.ll2011l) show that 
the NLTE effects strongly affect the less abundant species 
(i.e., Fei) which is over-ionized in RGs, in particular in 
metal-poor stars. This effect explain why the LTE ion- 
ization equilibrium between neutral and singly-ionized 
iron is destroyed and why we do not force [Fei/H] to be 
equal to [Fen/H] by decreasing the surfa ce gravity. Note 
that f or the five stars in common with Shetrone et al.l 
(2003) we have on average larger logg values by 0.2 dex, 
but they decreased their photometric estimates of logy 
by —0.29 dex. The decrease in the log <? values artificially 
forces the ionization equilibrium and causes an increase 
in their Fei abundance determination. This approach 
partially overcome the NLTE effects by changing the sur- 
face gravity. The more abundant species (i.e., Fen) are 
supposed to be free from NLTE effects, at least on the 
ionization equilibrium. 

5.4. The impact of NLTE effects 

To further constrain this effect, we directly compared 
Fei and Fen abundances. Data plotted in panels a), b) 
and c) of Fig. 6 indicate that the difference appear to in- 
crease when moving from metal-rich to metal-poor stars. 
We performed a linear fit (panel a, red dashed line, the 
zero-point a and the slope $ are labeled) over the 27 stars 
for which we have both Fe i and Fe n measurements. We 
found that the two metallicity distributions differ at the 
^97% confidence level. To further constrain this effect, 
we performed the same analysis, but using stars with 
abundances based on at least two Fen lines (panel b). 
The sample decreased to 19 stars, but the two distribu- 
tions still differ at the ^87% confidence level. The out- 
come is the same if we use stars with abundances based 

35 The star Car40 is slightly fainter than the selected stars 
(V=18.56), but the iron abundance is based on a single line. 



on at least three Fen lines (panel c), 15 stars, ~75% 
confidence level. 

On the basis of this empirical evidence, we corrected 
the Fei abundances to account for the NLTE effects by 
using the linear fit between Fei and Fen abundances 
plotted as a red dashed line in panel a) of Fig. 6. The new 
Fe i metallicity distributions accounting for the NLTE ef- 
fects are plotted as red dashed lines in panels a), c) and e) 
of Fig. 5. The weighted mean Fe i metallicity increases by 
~0.1 dex, while the weighted standard deviations attain 
similar values. The use of linear fits based on more accu- 
rate Fen abundances for either 19 (panel b) or 15 (panel 
c) stars gives similar changes in the metallicity distribu- 
tion J^l. These results indicate that the weighted mean 
metallicity of the Carina RGs, corrected for NLTE ef- 
fects, ranges from /x(NLTE) = -1.68 (15 stars) to -1.80 
(44 stars), while the weighted standard deviation ranges 
from 0.18 to 0.24 dex. The extreme range in iron abun- 
dance covered by the stars with the most accurate mea- 
surements is ~1 dex. 

In order to constrain the NLTE effects we performed 
NLTE computations for two metal abundances ([Fe/H] = 
—2.0, —1.5) using the Fei/n atom model provided by 
iCollet et alj ((2005). Note that to constrain the slope 
of the NLTE effects between Fei and Fen lines, the 
two marcs model atmospheres used in NLTE compu- 
tations have the same effective temperature (4500 K) 
and the same surface gravity (1.0 dex). The calcula- 
tio ns were performed follow ing th e approach adopte d 
by IThevenin fe Idiartl (fl~999h and bv ICollet et all ((2005T) . 
Moreover, we did not use inelastic collis ions for super- 
giants as suggested bv lMerle et alJ (|2011l ). so our errors 
on NLTE abundances have to be considered as upper 
limits. We found significant NLTE effects in both Fei 
and Fe n lines. Moreover, we found a significant rela- 
tive difference between [Fei/H] and [Fen/H] theoretical 
abundances. A glance at the data plotted in Fig. 7 dis- 
plays that the predicted slope agrees quite well with the 
observed one. In particular, the NLTE computations 
show stronger effects on Fei than on Fen abundances 
(~ 0.1 dex). However, more detailed calculations ac- 
counting for changes in effective temperature and grav- 
ity are required before we can reach firm conclusions con- 
cerning the impact of NLTE effects on Fe i and Fe n abun- 
dances. The current findings do not de monstrate but fur- 
ther s upport the results obtained by iBergcman n et all 
(2011) concerning the occurrence of NLTE effects in 
metal-poor RG atmospheres. Finally, we note that our 
preliminary calculations also show non-negligeable NLTE 
effects on the Fen lines. This finding, once confirmed by 
more detailed and independent calculations, implies that 
the true mean [Fen/H] abundance of Carina stars should 
be increased by at least 0.1 dex. 

6. FINAL REMARKS 

The current findings soundly support previous esti- 
mates of Carina's mean metallicity based on the differ- 
ence in color between the Red Clump (RC) stars and 
the horizontal-branch (HB) at mean color of RR Lyrae 

36 The weighted mean metallicities and the weighted standard 
deviations based on the linear fit plotted in panel c) of Fig. 6 are: 

(NLTE) = -1.75 ± 0.02, ct(NLTE) = 0.23 dex [44 stars], 
/i(NLTE) = -1.71 ± 0.03, ct(NLTE) = 0.24 dex [20 stars], 

(NLTE) = -1.64 ± 0.04, ct(NLTE) = 0.17 dex [15 stars]. 
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[Fe/H](A(f? — = -1.70 ±0.19 (IBono et al.ll2010h . 

The same outcome applies to the mean metallicity esti- 
mated usi ng stellar isochrones [Fe/H] ~ —1.79 (BaSTI 

data base. IStetson et aT] 120111 ). [Fe/H] 1.79 ± 0.35 

(jLianou et al.1120111 ). High- resolution spectra for 15 stars 
support, t ogether with s i milar results available in the 
literature (Lcmasl e et al.l I2012D . a significant decrease 
in the range in iron abundances when compared with 
similar measurements based on CaT (~ 1 vs 2-3 dex). 
The same conclusion applies to the comparison with the 
metallicity distribution predicted by chemical evolution 
models ba sed on star formation h istories available in the 
literature ([Lanfranchi et al]|2006l) . 

Our sample is more than a factor of four larger than 
any previous spectroscopic investigation based on high- 
resolution spectra. However, the current data do not 
allow us to determine whether the spread is either atmo- 
spheric, i.e., caused by a difference in the mean metallic- 
ity between the old and the intermediate-age population, 
or by measurement errors. To assess whether the differ- 
ent stellar populations are also characterized by different 
mean metallicities, new spectra with high S/N down to 
RC (intermediate-mass) and to red HB (low-mass) stars 
are required. 

We determined the mean iron abundance for a sample 



of 44 Carina RG stars. We corrected for NLTE effects 
on the ionization equilibrium and we found [Fe/H] = 
— 1.80 and a — 0.24 dex. If this estimate of Carina's 
mean metallicity is supported by future, the position 
of this dSph in the Metallicity-Luminosity diagram will 
be between one and two a more metal-rich than ex- 
pected according to the empirical relation followed by 
dSph, dE/dSO and giant early type galaxies (|Mateoll20"08t 
Chilingari an et al.ll2011[ ). If this turns out to be the case, 
it might open new issues concerning the interplay be- 
tween chemical evolution and enrichment, star formation 
history and stellar evolution in gas-poor ste llar systems 
(jWalker et alJ l2009: R evaz fe Jab lonka 20i3). 
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APPENDIX 

The spectroscopic target Car3 that we have in c ommon with the sa mple co llected bv | Shetrone et all (|2003l Car3) 
deserves a more detailed discussion. According to IBono et ail ()2010fl and to IStetson et al.l ()201lD this star has the 
following photometry: U = 20.700 ± 0.033, B = 19.172 ± 0.001, V = 17.675 ± 0.003 and J = 16.102 ± 0.003 mag. The 
position of the star in the U, B~I CMD suggests that it might not be a Carina member, but a foreground field star. 
The V, B-I CMD and the B-V, V-I color-color diagram support the same evidence, but not as strongly. The star 
has five fainter companions within 5", all of them with AV > 4.8 mag, and a distance > 2.4". The star is located 2.5' 
N-NE of the galaxy center. The neighborhood stars are sufficiently faint and distant that the luminosity contamination 
on Car3, with a typical seeing of 1", is at most of the order of one part in 10 4 . We also checked individual photometric 
measurements and we found marginal evidence of variability. 

The radial velocity of Car3 is 231.0 ± 0.8 km s" 1 , i.e. v ery similar to the mean radial velocity of candidate Carina 
stars (RV mcan = 220.9 ± 0.1 km s" 1 , iFabr izio et all 1201 ID . To constrain whether Car3 is a giant or a dwarf we also 
forced the balance between Fe i and Fe n lines and we found that the new estimate of the surface gravity is within 
0.1 dex from the estimate based on the photometry (logg = 0.45±0.19 dex). This finding further support the evidence 
that Car3 is a t ruly red giant. 

According to iShetrone et al.l (|2003l ) Car3 is underabundant in a-elements, and in particular in Ca. To further 
constrain this evidence we selected three Carina RGs (Car4, Car37, Car50) with similar surface gravities (logg = 
0.53 to 0.67 dex) and effective temperatures (T e s = 4220 to 4290 K). We performed a detailed Ca abundance analysis 
by using two diffe rent multiplets, n amely A = 6102, 6162, 6162 A and A = 6166, 6169.0, 6169.5 A. We also accounted 
for NLTE effects (jMerle et al.ll201lD and we found that Car3 is on average a factor of two more Ca poor that the other 
three Carina RGs. We also visually inspect the entire spectrum and we did not find evidence of mo lecular bands, thus 
suppo rting the absence of this star in the list of Carina carbon stars provided bv lMould et al.l (|1982l) and bv lKoch et al.l 
(2008). Interestingly enough, we found that Car3 shows a strong H a line with both a blue and a red P-Cygni profile. 
Thus, suggesting that this star might have a complex atmosphere, probably affected by chromospheric activity. Car3 
is an object that deserves further photometric and spectroscopic investigations. 
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TABLE 2 

Atomic Data and Individual Iron Abundances 1 











Car2 




Car3 




Car4 




CarlO 




A (A) 


Elcm. 


EP (eV) 


logs/ 


[Fe/H] 


EW a 


[Fe/H] 


EW a 


[Fe/H] 


EW a 


[Fe/H] 


EW a 


A HO A T7C\ 


26.0 


2.279 


—2.241 


















a n cr n 1 1 r\ 

4y5u.nu 


26.0 


3.417 


— 1.670 


















a r\To "i nn 

4y73.10U 


26.0 


3.960 
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on 1 r\ 1 1 

— z.zU ± U.14 


88 










c n /i a 01 n 

5U44.21U 


26.0 


2.851 


—2.038 






— 2.15 ± U.13 


93 






— 2.10 ± U.12 


51 


5U54.64U 


26.0 


3.640 


— 1.921 


















5090.770 


26.0 


4.256 


-0.400 


















5197.580 


26.1 


3.230 


-2.348 


-1.55 ± 0.30 


105 


-1.85 ±0.38 


117 


-1.49 ±0.42 


105 


-2.00 ±0.21 


73 


5215.180 


26.0 


3.266 


-0.871 


















5217.390 


26.0 


3.211 


-1.070 


















5234.620 


26.1 


3.221 


-2.279 






-1.95 ±0.34 


83 










5242.490 


26.0 


3.634 


-0.967 






-2.05 ±0.17 


98 










5264.810 


26.1 


3.230 


-3.133 


















5284.110 


26.1 


2.891 


-3.195 










-1.55 ±0.20 


65 






5288.520 


26.0 


3.694 


-1.508 


-1.65 ± 0.13 


52 






-1.70 ±0.24 


56 






5307.360 


26.0 


1.608 


-2.987 


















5316.620 


26.1 


3.153 


-2.014 






-2.05 ±0.20 


156 











This table is available entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content. 

a Equivalent widths are in mA. 



TABLE 3 

Impact of uncertainties in effective temperature, surface 
gravity and microturbulence velocity on iron abundances for 
the representative star car23. 





AT cff 


[K] 


A log g 


log(cms 2 )] 


A£ [km 


s- 1 ] 




Ion 


-100 


±100 


-0.3 


+0.3 


-0.5 


+0.5 


<o-> a 


Fcl 


+0.10 


-0.11 


-0.01 


-0.00 


-0.10 


-0.01 


0.11 


Fell 


±0.02 


±0.10 


-0.02 


+0.13 


-0.09 


-0.04 


0.12 



Weighted standard deviation. 
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Fig. 1. — Top: V, B-V CMD of the spectroscopic targets. The black dots mark the UVES targets with a RV between 180 and 260 
km s — 1 (Carina candidate members), while the small orange dots show targets with RV outside the previous range (field star s). The two 
large o rang e dots refer to the two carbon stars. The blue squares and red stars mark the ta rgets in common with the sam ple of Ko ch et al.l 
(2008) and Shctronc et al. (2003). The colored lines show three scaled Solar isochrones (Pictrinfcrni et al. 2004, 2006) at the same age 
(t = 12 Gyr) with [Fe/H] = -1.50 in red, [Fe /H] = -1.79 in black and [Fe/H] = -2.27 in blue. The background gray dots are from the 
Carina photometric catalog (Bono et al. 2010). Bottom: Same as the top, but the targets are plotted in the effective temperature vs surface 
gravity plane. The error bars display typical uncertainty on both effective temperature and surface gravity. 
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Fig. 2. — Detailed iron abundances as a function of the Excitation Potential (EP) for five stars covering the range in visual magnitude of 
the entire sample. From top to bottom the panels show for each star the line-by-line analysis for the Fe I abundance (left panels) and for 
the Fell abundance (right panels). For each ionization stage, the mean weighted value and the statistical error are given and plotted as a 
dotted line. Numbers in square brackets are the number of iron lines used to measure the abundance. 
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Fig. 3. — Panels a): Difference A[Fe/H] = [Fc/H] our — [Fe/H] author in Fel (left) and in Fell (right) stellar abundances between our 
measurements and similar measurements by Shetrone ct al. (2003). The weighted mean fi (thick lines) and the weighted standard deviation 
a are labeled together with the 1 a interval (dotted lines). The error bars display i ndividual uncertai nties of the two sets of abundances. 
Panels b): Same as panels a), but the difference i s with similar measur ements from Koch ct al. (2008). Panels c): Same as panels a), but 
the difference is with similar measurements from Lemasle et al. (2012). Panels d): Same as panels a), but the difference is with similar 
measurements from Venn et al. (2012). 
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Fig. 4. — Difference in iron abundance (as in Fig. 3) for the 25 Carina stars in common with Koch ct al. (2006). They estimated stellar 
iron abundances from the CaT measurements based on medium-resolution spectra in the ZW84 (left) and in the CG97 (right) metallicity 
scale. 
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Fig. 5. — Panels a)— b): Metallicity distributions of Carina stars based on Fel (left) and Fell (right) lines. The weighted means and the 
weighted standard deviations are labeled, together with the sample size. The red dashed line shows the Fe I metallicity distribution of 
the same stars, but accounting for NLTE effects. The correction of individual abundances is based on the linear fit between Fel and Fell 
abundances given in panel a) of Fig. 6. The new /x(NLTE) and ct(NLTE) values are labeled in red. Panels c)-d): Same as the top, but for 
stars with more than twelve Fel or with at least two Fell measurements. Panels e)— f): Same as the top, but for stars with at least three 
Fell measurements. 
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Fig. 6. — Panel a): Comparison between stellar iron abundances based on neutral and on singly-ionized lines. The red dashed line shows 
the linear fit of the data, used to correct the Fel metallicity distribution for NLTE effects. The zero-point (a) and the slope of the 
linear fit are labeled. Panels b) and c): Same as panel a), but for stars with at least two and at least three Fell measurements, respectively. 




Fig. 7. — NLTE effects on ionization equilibrium. The symbols (black dots) are the same as in panel a) of Fig. 6. The blue diamonds are 
the abundances given by NLTE computations of Fel/II for two representative models of Carina RGs. We adopted an effective temperature 
of 4500 K, a surface gravity of 1.0 dex and two iron abundances [Fe/H] = —2.0, —1.5. The dashed red line shows the linear regression on 
the data (the same as in panel a of Fig. 6), while the solid blue line shows the linear fit on the two NLTE models. The zero-point (a) and 
the slope (/3) of the linear fits are labeled. 



